Electronic phase diagrams of carriers in self-assembled InAs/GaAs quantum dots: 
violation of Hund's rule and the Aufbau principle for holes 
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We study the orbital and spin configurations of up to six electrons or holes charged into self- 
assembled InAs/GaAs quantum dots via single-particle pseudopotential and many-particle config- 
uration interaction method. We find that while the charging of electrons follows both Hund's rule 
and the Aufbau principle, the charging of holes follows a non-trivial charging pattern which violates 
both the Aufbau principle and Hund's rule, and is robust against the details of the quantum dot size. 
The predicted hole charging sequence offers a new interpretation of recent charging experiments. 

PACS numbers: 73.21. La, 73.23.Hk 73.63.Kv 
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The remarkable combination of three-dimensional spa- 
tial confinement in quantum-dots with the ability to in- 
tegrate them into carrier-transporting device structures 
enables storage and retrieval of electrons 0, S 0, @ , 
to the benefit of future quantum-computing [3, |8j , mem- 
ory 0,01, and single-photon |0 applications. Unlike 
real atoms, where large Coulomb repulsion energies J ~ 
10 eV limit the number of ionized species to just a few, 
semiconductor quantum dots can be loaded by as many 
as six Q to ten [1^ electrons in colloidal T2\ and self- 
assembled P, 13,01 dots, and up to hundreds of electrons 
in electrostatically confined dots Furthermore, 
one can measure for each ionization state the stable spin- 
configuration 111 II, 1MB, the energy to add another elec- 
tron 0, 0, H-Tsl I12II as well as the attendant spectro- 
scopic shifts with charging 0,^3- Like real atoms, the 
stable spin configuration observed in electrostatic dots 
HHH, having lateral dimensions of 500-1000 A, follow 
the rules of atomic physics; that is, the s, p, d, ... shells 
are occupied in successive order with no holes left be- 
hind (Aufbau principle) and with maximum spin (Hund's 
rule). Recently, it became possible to load and measure 
electrons 0, and holes 0, 0, into much smaller 
(~ 200 X 40 A), epitaxially grown self-assembled dots 
of InGaAs/GaAs. The analysis of the results, however, 
was based on the effective-mass approximation (EMA) 
theoretical framework [13], (which is designed for dots 
much larger than the exciton Bohr radius), leading to 
the conclusion that electronic configurations are atomic- 
like. Electronic structure calculations for self-assembled 
dots ^3 reveal that while for electrons the Coulomb en- 
ergy Jee ~ 20 meV is smaller than the level spacing Ae ~ 
50-70 meV, for holes the Coulomb repulsion J^h — 15-25 
meV is comparable to the level spacing Ae ~ 10-20 meV. 
This opens the possibility of observing for holes some 
stable, exotic spin configurations that defy the rules of 
atomic physics. We have applied a combination of an 
atomistic pseudopotential description (lij for the single- 
particle level structure, with a many-body configuration 
interaction (CI) description 0| of many-particle effects 
to both electron and hole loading into InGaAs/GaAs self- 



assembled quantum-dots. We calculate the generalized 
electronic phase diagram of the system showing which 
many-particle configurations are energetically stable for 
various p-p and p-d splitting of the single-particle levels. 
We find that while electron loading follows both the Auf- 
bau principle and Hund's rule, hole loading gives rise to 
stable but unusual spin configurations. While these cal- 
culated configurations agree with recent measurements 
[1^ lisl , they differ from their interpretation Hsf. which 
assumes 2 dimensional (2D) parabolic models [19j that 
have been employed extensively and successfully to ana- 
lyze large electrostatically confined dots. ^ H, 01 - The 
reason for the failure of the simpler interpretation is that 
parabolic models ignore the inter-band and inter-valley 
coupling existing in a real self-assembled quantum dot. 

The "charging energy" fi{N) is the energy needed to 
add a carrier to the dot that is already charged by — 1 
carriers, fj,{N) = E{N) - E{N - 1), where E{N) is the 
correlated many-body total energy of the ground state of 
the A^-particle dot. The "addition energy" A{N - 1,N) 
(analogous to the difference between ionization potential 
and electron affinity) indicates how much more energy 
is needed to add the A^th carrier compared to the en- 
ergy needed to add the A^ — 1th carrier: A(A^ — 1, A^) — 
^l{N)-fi{N~l) = E{N)-2E{N-l) + E{N-2). In the 
Hartree Fock approximation, where the effect of correla- 
tions is neglected but the direct Coulomb and exchange 
interactions are retained, simple expressions can be de- 
rived for the addition energies. To do so, one needs to 
decide first what is the order of filling the single-particle 
s, p d, ■ ■ ■ levels. Assuming the filling sequence of the 
left (L) hand side of Fig. obeying Hund's rule and 
the Aufbau principle, the addition energies are given by: 
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(b) Reuter's model 
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(c) Pseudopotential-CI 
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FIG. 1: (Left) The ground state configurations (a) calculated from 2D parabolic model, (b) suggested by Reuter et al and 
(c) calculated from pseudopotential-CI method. (Right) The phase diagrams for A'^ = 4, 5, 6 electrons/holes (a) calculated from 
2D parabolic model, (b) for electrons using realistic Coulomb and exchange integrals, and (c) for holes using realistic Coulomb 



and exchange integrals. For A'^=6, 5di,d2 = ^vi,V2 assumed. 

self-assembled dots. The labels "e" and "h" inside the symbols denote electron and hole respectively, 
phases I, II, III, and IV are given in the text. 



The circles represent electrostatic dots, while lens represent 

The configurations of 



where Jij and Ki^j are, respectively, the Coulomb and 
exchange integrals between states i and j . 

To calculate these addition energies, one must first con- 
struct a single-particle Schrodinger equation model. In 
this step, one might need to account not only for quan- 
tum confinement, but also for electronic structure effects 
such as multi-band (light-hole, heavy-hole, conduction) 
coupling; inter-valley iV-X-L) coupling; spin-orbit cou- 
pling; and the effect of strain and chemical intermixing. 
It is then possible to compute all of the single-particle 
level spacings and integrals entering Eq. J^l, thus pre- 
dict the value of Ahf(-^ — 1, N). Alternatively, one can 
neglect explicitly electronic structure effects other than 
quantum confinement, and use instead a particle- in-a- 
parabolic-box model, widely used in this field [iH Il9| . 
In this 2D-EMA, the p levels are degenerate (epi=ep2), 
as are the d levels (cdi = = ^ds) ^-nd the splitting 
between the s and p levels (Es — Sp) and the splitting 
(£p — £d) between the p and d levels are both equal to 
the harmonic oscillator frequency uj. Furthermore, the 
assumed parabolicity assures that analytic formulas can 
be derived fo^' the Coulomb and exchange matrix el- 
ements that relate all integrals needed for the addition 
energies to the value of a single Jss, for instance in Eq. 

2Jsp — Jss — Ksp — Jss/ 4:- Thus, universal results 
can be derived for electrons and holes as shown in the 
right (R) hand side of Fig. for 7V=4, 5, 6. Since the 



restriction of the 2D-EMA model to degenerate shells 
(epj=ep2, edi — = ^da ) and to equidistant shells 
(cp — es = Cd — ep) might be rather stringent [20j. we allow 
in Fig. [TJi(R), Sp^^p^ = e^^ - ep^ and Sp^^di = Cdi - to 
vary, calculating for each {N; Sp-^^p^, '^^2,^1} the configura- 
tion which minimizes the total-energy. This gives a phase 
diagram as a function of the parameters Sp-^^p^ and Sp^.di 
in the unit of Jss as shown in Fig. ^(R) for particle num- 
ber 7V=4, 5 , 6. The 2D-EMA model yields for N^A two 
electronic phases: high-spin l4—{s^s^){p\){pl) and low- 
spin ll4={s'^ s^){p\p\). For = 5, we find three phases 
15= {s^s^){pl){pl){dl), Il5= {s^ si){plpi){pl) and III5- 
{s^ s^){plp\){d\). For A^ = 6, we find four phases, le^ 

isU^){pl)ipl)idl)idl), Il6= {sh^)iplpi)ipl){dl), IIl6 = 

ish^)iplpi)iplpi) andIV6= {sh^)iplpi)idldi). Tode- 
cide which of these phases is a ground state, we need to 
know in Fig. ^(R) the value of 5pi^p2/Jss and Sp2.di/ Jss- 
For electrons in self-assembled dots, the single-particle 
energy spacing is usually more than twice the Coulomb 
energy, so, Sp2,di > 2 Jss Q, El- For holes, i5 „ 2. rf i = 
1.17 Jss was determined from recent experiments [l^ [iSj 
and (5pi,p2 = is assumed in 2D-EMA model. This places 
in Fig. ^(R) for both electrons and holes, phases I4, II5, 
Illg as ground states for A^=4, 5, 6, respectively. The 
ground state configurations of the 2D-EMA model are 
collected in Fig. [Tli(L), for N=l - 6. 

For electrons, the ground states of 2D-EMA model 



3 



TABLE I: Hole addition energies of self-assembled InAs/GaAs quantum dots in meV. The experimental values are from 
Ref. |_15j at zero magnetic field. The "Ground State" values correspond to the low-spin configurations as given in Fig. ^(L) 
and the "Excited State" values to the high-spin configurations assumed in Ref. s^nd given in Fig. EKL). The results of 
"Pseudopotential-t-CI" calculations correspond to the configurations from Fig. 0;(L). 
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are corroborated by atomistic pseudopotential calcula- 
tions [Fig^(R)], where we use the Coulomb integrals 
obtained from atomistic wavefunctions for electrons in a 
lens shape InAs/GaAs dot with 25 nm base and 2.5 nm 
height. This shape is realistic, according to experimental 
findings |2ll |. and predicts a fundamental photolumines- 
cence line very close to the one observed in the charging 
experiment jl5l | at around 1 eV. Overall, the compari- 
son between Figs. 2a and Figs. 2b shows that while the 
phase boundaries can change significantly when realistic 
wavefunctions are assumed instead of 2D-EMA values, 
the ground state symmetries for N=5 and 6 electrons in 
self- assembled dots remain unchanged and are far from 
other competing phases. 

The foregoing analysis of loading of electrons ^ H, IS 
has been simplified by the fact that the single-band ef- 
fective mass model is not a drastic approximation given 
that in direct-gap zinc-blende materials electrons derive 
from a nondegenerate, spin-orbit-free Fic band which 
is energetically isolated from other states. The analy- 
sis of loading of holes, however, does not benefit from 
the same simplification, as holes derive from a mix of 
heavy- and light-hole states, invalidating ji^ the classi- 
fication of hole states as pure s or p or d levels and as 
pure heavy- hole or light-hole states Furthermore, 
as shown in Fig. ^(R), unlike electrons, the ground 
states of holes are close to competing phases in the 2D- 
EMA phase diagrams, which require a more careful treat- 
ment. Nevertheless, the 2D-EMA model is still attrac- 
tive in its simplicity and Renter et al. [Tsf used it to 
analyze their hole charging results. The value of Jss is 
directly accessible from experiments since it is well ap- 
proximated by A(l,2). The only remaining parameter 
in the 2D-EMA model is the single-particle energy split- 
ting (eg — ep)—{tp — td)=^ which can be extracted from 
measuring A(2, 3) = (ep - e^) + Renter et al [l3| 

thus determined Jss— 23.9 meV and Sp — Ss = 28 meV. 
Since experimentally five addition energies are available 
and only two were used in the fit above, the problem is 
overdetermined and it is possible to assess how well the 



model reproduces the remaining three experimental data 
points. Assuming 5pi^p2—0 and Sp2,di=^d — <^p ~ip~^s—28 
meV, yields the above mentioned Sp2,di~^-^7 Jss, leading 
to the phase diagram of Fig. ^^^(R), with ground state 
phases: I4, II5 and Ille- If one calculates the addition en- 
ergies of Eq.(P), using these ground-state configurations 
one gets the values indicated by "2D-EMA model ground 
state" in Table However, this hole addition sequence 
contradicts the magnetic field data of Ref. Iisl which 
show that the hole c?-levels are occupied before p-levels 
are filled completely (non- Aufbau}_|l5j . To explain their 
magnetic field data. Renter et al. [l^ assumed a hole fiU- 
mg sequence [Fig.[Tl3(L)] that relies on an ad hoc excited 
hole state instead of the the ground states predicted by 
the 2D-EMA model [Fig.[Tli(L)]. In Tabled we compare 
the ensuing addition energies of both the ground state 
and excited state configurations with experiments. We 
find that the addition energies given by both scenarios of 
Fig. ni(L) and^(L) show significant discrepancies from 
the experimental values, with about 25 -50 % error. 

For holes, the discrepancy between the harmonic os- 
cillator results Fig. ^^^(L) (Aufbau configuration) and 
the magnetic field data of Fig. ^^(L) (non-Aufbau con- 
figuration) is difficult to reconcile within the 2D-EMA 
model, since even relative large variations of the single- 
particle energy parameters <5pi,p2 and 6p2^di do not lead 
to a ground-state configuration change from Fignji(L) 
to Fign]3(L). Our atomistic pseudopotential, plus CI 
calculations show different ground-state configurations 
for holes than the above two models. As shown in 
Fig^(R), for A^=4, and iV=5, using an atomistic de- 
scription the topology of the phase-diagrams are the 
same as in the 2D-EMA model, but the boundaries are 
shifted. As a result, for A^=4, the hole ground state 
is now phase ll4={s'^ s^){plp{), not ii—is^ s-^){pl){pl) 
, and for iV=5, the ground state is phase IIl5 = 
{sls^){plpi)idl), not 115= ish^)iplpi)ipl). For A^=6, 
the topology of the phase digram changed completely: 
phase l6={s^ s^){pl){pl){dl){dl) disappeared, and the 
ground state is now IV6= {s'^ s^){plp{){dld{). The 
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TABLE 11: First hole Coulomb energy Jss and single-particle 
energy level spacings in meV, from atomistic pseudopoten- 
tial calculations for six different self-assembled len-shaped 
InAs/GaAs quantum dots. 
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ground-state configurations are listed in Fig^(L) for 
N=l - 6. Using these new ground states, Table com- 
pares the experimental addition energies and the cal- 
culated results for six different InAs/GaAs lens shaped 
quantum dots of different bases and heights. Very good 
agreement is obtained for the InAs dot with 20 nm base 
and 2.5 nm height, with differences in the addition ener- 
gies of less than 16%, compared with almost 50% error 
in 2D-EMA model (despite the fact that two of the ad- 
dition energies were fitted). The parameters Spi^p2 and 
Sp2,di calculated for different dots are given in Table UTI 
and as shown in Fig. ^^(R), lie close to the center of the 
predicted phases II4, III5 and IVe. This indicates the 
stability of our numerical results against possible varia- 
tions of 6pi,p2 and Sp2.di due to shape anisotropy or alloy 
effects. The addition energies we calculate from config- 
uration interaction are different by about 1-3 meV from 
those we calculate in Hartree-Fock, although they give 
the same ground state. 

Our predicted charging pattern [Fig. ^^(L)] shows that 
the level filling by holes does not follow the Aufbau prin- 
ciple nor the Hund's rules: d levels get filled before the 
second p level, despite the fact that the d level is en- 
ergetically more than 3 meV above the second p level. 
The non-trivial hole filling pattern is due to two reasons. 
First, the large p-level splitting leads to the p2 level being 
energetically close to the di level, i.e., Sp^^di is small. In 
Table we list the hole single-particle energy spacings 
and the first hole Coulomb integrals Jss for different dots. 
We see that Sp^.p^ ~ (0.3 - 0.5) -Jss and Sp^.di ^ (0.2 - 
0.3) -Jss, which differ significantly from the assumption 
of 2D-EMA where, Spi,p2 ~ 0, and 5p2,di ^ Jss- Second, 
the Coulomb repulsion between the pi and the d level is 
lower than that between the two p levels, therefore the 
Coulomb energy can overcome the single-particle energy 
spacing Sp^^di, leading to the non- Aufbau charging pat- 
tern. 

An important feature of the present theory is not only 
its compatibility with the zero field experimental results 
but also with the magnetic field dependence obtained in 
Ref. 0. The experimental results from Ref. show 
that the spin Zeeman effect is significantly smaller than 
the orbital Zeeman and in the interpretation of the ex- 



perimental results only the latter has to be taken into 
account. Within this approximation our theoretical re- 
sults reproduce the experimental data for iV = 1 - 6 holes 
given in Fig. 2 of Ref. I15L since our results have the same 
total orbital angular momenta. The difference in the to- 
tal spin for configurations assumed in Ref. and in 
the present work is significant. We suggest to identify 
the spin configurations by measuring optically the fine- 
structure of the N h + e ^ {N ~ 1) h transition in large 
magnetic field. While our suggested configurations for 
the TV = 4, 6 holes have closed-shell fillings and only one 
single line with no fine-structure should be observed, the 
high-spin configuration suggested by Renter et al. [Tsf 
should show a rich fine-structure that should become ap- 
parent, especially at high magnetic fields. 

In conclusion, we analyzed the many-particle config- 
urations for electrons and holes in quantum dots in the 
form of phase diagrams. From these diagrams, we pre- 
dict that the hole charging sequence presents surprising 
configurations (not expected from effective mass calcula- 
tions) that violates the Aufbau principle as well as the 
Hund's rule. Our results are in good agreement with re- 
cent experimental findings and provide a novel way to 
study the charging of carriers in quantum dots. 

This work was supported by the US Department of 
Energy, Office of Science, Basic Energy Science. 
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